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ABSTRACT 


A mathematical model is developed to predict the dynamics of the proposed 
Spacecraft Control Laboratory Experiment during the station keeping phase. The 
Shuttle and reflector are assumed to be rigid, while the mast connecting the 
Shuttle to the reflector is assumed to be flexible with elastic deformations 
small as compared with its length. It is seen that in the presence of gravity- 
gradient torques, the system assumes a new equilibrium position primarily due 
to the offset in the mast attachment point to the refector from the reflector’s 
mass center. Control is assumed to be provided through the Shuttle’s three 
torquers and through six actuators located by pairs at two points on the mast 
and at the reflector mass center. Numerical results confirm the robustness 
of an LQR derived control strategy during station keeping with maximum control 
efforts significantly below saturation levels. The linear regulator theory 
is also used to derive control laws for the linearized model of the rigidized 
SCOLE configuration where the mast flexibility is. not included. It is seen 
that this same type of control strategy can be applied for the rapid single 
axis slewing of the SCOLE through amplitudes as large as 20 deg. These re^:- 
sults provide a definite trade-off between the slightly larger slewing times 
with the considerable reduction in over-all control effort as compared with 
the results of the two point boundary value problem (TPBVP) application of 
Pontryagin's Maximum Principle. In connection with the TPBVP, a useful numer- 
ical solution procedure for minimum time, 2-D and 3-D attitude maneuvers of 
a rigid spacecraft is developed and applied to the SCOLE system. Both singu- 
lar and nonsingluar cases can be handled. The minimum time is determined 
by sequentially shortening the slewing time. It is shown that one of the 


four initial costates associated with the quaternions can be arbitrarily 
selected without affecting the optimal controls, thus resulting in a simpli- 
fication of the computation. Finally, a stability criterion for a controller 
which operates in the continuous time domain but depends on discretized input 
data is developed and relates the maximum tolerable discretization step size to 
the damping ratio and undamped frequency of any mode in the system model. Com- 
putational requirements for the estimator and controller are evaluated based 
on a type 80387 microprocessor, and assuming that the number of actuators 
and sensors are a certain fraction of the number of state components. 
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INTRODUCTION 


The present grant, NSG-1A14, Supplement 10, represents an extension to 
the research effort initiated and accomplished in previous grant years (May 
1977 - March 1987) and reported in Refs. 1-14*. The attitude and shape con- 
trol of very large inherently flexible proposed future spacecraft systems is 
being investigated. Possible future applications of such large spacecraft 
systems (LSS) include: large scale multi-beam antenna communication systems; 
Earth observation and resource sensing systems; orbitally based electronic 
mail transmission; as platforms for orbital based telescope systems; and as 
in-orbit test models designed to compare the performance of flexible LSS 
systems with that predicted based on computer simulations and/or scale model 
Earth-based laboratory experiments. In recent years the grant research has 
focused on the orbital model of the Spacecraft Control Laboratory Experiment 
(SCOLE) first proposed by Taylor and Balakrishnatf^ in 1983. 

The present report is divided into six chapters. Chapter II is based on 
a paper presented at the 1988 AIAA/AAS Astrod 5 nnamics Conference which describes 
the development of a mathematical model to predict the dynamics and control 
of the SCOLE configuration in-orbit during station keeping. The flexibility 
of the mast which connects the reflector to the Shuttle is included in the 
model. In accordance with the SCOLE design challenge^^ the attachment point 
between the mast and the reflector is considered to be offset in two cartesian 
directions from the mass center of the reflector. The linear regulator theory 
is used to derive both orientation and mast vibration suppression control laws. 


* References cited in this report are listed separately at the end of each 
chapter. 
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In the following chapter (Chapter III) the same linear regulator theory 

techniques are applied to develop control laws of a rigid model of the SCOLE 

which could be used for rapid single-axis slewing through as much as 20 deg. 

These control laws could provide an attractive alternative to the slewing 
14 

strategies based on the two point boundary value problem associated with 
Pontryagin's Maximum Principle where a trade-off could be made between the 
noticeable savings in control effort and the increased slewing times associ- 
ated with the LQR applications. Chapter III is based on a paper accepted 
for presentation at the 39th International Astronautical Congress in 
October 1988. 

Chapter IV represents a completely revised version of a paper presented 
at the AIAA 26th Aerospace Sciences Meeting and describes a numerical approach 
for solving general three dimensional rigid spacecraft minimum time attitude 
maneuvers. In place of the total slewing time an integral of a quadratic 
function of the controls is used as the cost function and allows for the 
treatment of both the singular and nonsingular problems in a unified way. 

The resulting two point boundary value problem is developed by applying the 
Maximum Principle to the system and solved by using a quasilinearization 
algorithm. 

In the following chapter (Chapter V) aspects of the computational 
requirements for the implementation of the controller are discussed. An 
improved stability criterion for a controller, designed to function in the 
continuous time domain, but which receives discretized observational inputs, 
is developed. An expression for the maximum tolerable discretization step 
size is derived in terms of the undamped frequency and damping ratio of 
any mode in the continuous-time system model. In the second part of Chapter V 
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the (on-board) computational requirements for the estimator and controller 
are evaluated based on a type 80387 microprocessor, and assuming that the 
number of actuators and sensors are a certain fraction of the number of 
state components. 

Finally, Chapter VI describes the main general conclusions together 
with future recommendations. The effort described here is being continued 
during the 1988-89 grant year in accordance with our most recent proposal. 
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II. 


THE DYNAMICS AND CONTROL OF THE ORBITING 
SPACECRAFT CONTROL LABORATORY EXPERIMENT 
(SCOLE) DURING STATION KEEPING 


ORIGINAL PAGE IS 
OF POOR QUALITY 


Abstract 


A mathematical model is developed to predict 
Che dynamics of the proposed orbiting Spacecraft 
Control Laboratory Experiment during the station 
keeping phase. The Shuttle as well as the reflec- 
tor are assumed to be rigid, the mast is flexible 
and is assumed to undergo elastic displacements 
very small as compared with its length. The equa- 
tions of motion are derived using a Newton-Euler 
formulation. The model Includes the effects of 
gravity, flexibility, and orbital dynamics. The 
control Is assumed to be provided to the system 
through the Shuttle's three torquers, and through 
six actuators located by pairs at two points on 
the mast and at the mass center of the reflector. 
At each of the locations, an actuator acts par- 
allel to the roll axis while the other one acts 
parallel to the pitch axis. It Is seen that. In 
the presence of gravity-gradient torques in the 
system dynamics, the system assumes a new equili- 
brium position about which the equations must be 
linearized, primarily due to the offset in the 
mast attachment point to the reflector. The 
linear regulator theory Is used to derive control 
laws for the linear model of the SCOLE Including 
the first four flexible modes. Numerical results 
confirm the robustness of this control strategy 
for station keeping with maximum control efforts 
significantly below saturation levels. 


I. Introduction 


The problem of maneuvering a flexible space- 
craft while suppressing the Induced vibrations Is 
becoming Increasingly Important. NASA Is Involved 
In studies which are concerned with the. control 
of flexible bodies carried by a Shuttle in an 
Earth orbit. Similar experiments are being con- 
ducted In Earth-based laboratories. It Is then 
desirable to derive a formulation which can accommo- 
date both types of experiments. 

NASA is currently involved In at least two 
experimental programs to test techniques derived 
for active control of flexible space structures. 

In several versions of a recent paper, SCOLE^ 
(Spacecraft Control Laboratory Experiment), Lawrence 
W. Taylor, Jr. and A.V. Balakrishnan have described 
the first which Is ground based. It Is a labora- 
tory experiment based on a model of the Shuttle 
connected to a flexible beam with a reflecting 
grillage mounted at the end of the beam (Fig. 1). 


As a part of the design challenge, the authors 
stressed the need to directly compare competing 
control desig;n techniques and discussed the feasi- 
bility of such a direct comparison. Concern would 
be given to modeling order reduction, fault manage- 
ment, stability, and dynamic systems. The second 
experimental program Is known as Control of Flex- 
ible Spacecraft (COFS)^ and consists of experi- 
ments designed to control flexible bodies carried 
by a Shuttle In an Earth orbit. Because of the 
cost and risks Involved In testing control techni- 
ques In space, COFS includes laboratory simulations 
of similar experiments which will precede the space 
test. Therefore, In assuring the success of both 
SCOLE and COFS, mathematical modeling and conq>uter 
simulation are required. 


To accurately model and simulate flexible 
spacecraft , one needs a thorough knowledge of its 
structural behavior. In a paper^, subsequent to 
the design challenge, the modal shapes and fre- 
quencies for the SCOLE system were derived. In 
reference 3, the SCOLE system Is assumed to be 
described by partial differential equations In 
which the variables separate. 


Undertaken In this study Is the modeling of 
the three dimensional dynamics of the SCOLE con- 
figuration based on the Eulerlan technique. This 
consists in isolating an elemental mass of the sys- 
tem In Its deformed state and deriving its angular 
momentum taken at the mass center of the Orbiter. 
The position vector extending from the origin of 
the coordinate syst€»n to the elemental mass of the 
mast or the reflector accounts for the elastic dis- 
placements. The expressions for these displace- 
ments are derived from the mode shape functions 
generated during the three dimensional structural 
analysis of the system. 

The equations obtained for the elemental mass 
of the components of the system are Integrated over 
the mass of the entire system to yield its angular 
momentum about the mass center of the Orbiter. The 
derivative of the system angular momentum with re?* 
spect to time Is equated to the gravity^radlent 
(and other external) torques on the system about 
the same point. Such a vectorial equation, when 
projected along the three axes of rotation, yields 
the system rotational equations of motion. These 
rotational equations of motion are then linearized 
to yield a model which provides the basis for the 
control law synthesis developed in this study. 


II. Angular Momentum of the SCOLE System 

A. Angular Momentum of the Shuttle about Its Mass 
Center, G 

The angular momentum of the Shuttle, taken as 
a rigid body, about its center of mass, G is 


2.1 


”s/G “ ^S/G ^S/Rq 

m 

where is the inertia tensor of the Shuttle and 

^S/Rq Shuttle’s inertial angular velocity. 

B. Angular Momentum of the Beam About G 

Consider an element of mass, dm, of the beam 
located at some point, P, such that (Fig. 1) 


^ M r r 

^/c * Y, {[fl^cos (wt+a)f^-<D sin («t+Y)£2 

+ [<» sin cos (wt+Y)f 2 -nyL^] j 

+ cos (u)t+a)f^+fty cosCwt+Y)!^] k} (6) 

C. Angular Momentum of the Reflector About, G 


GP - rQ + q (2) 

A 

where Tq ■ -zk is the position vector of P in the 
underformed state; q (z,t) * u(z,t) i + v(z,t)j 
in which, u and v are the x and y components of 
the mode shape vector, respectively. 

The angular momentum of dm about G, is 

given by; 

where 

r * -zk + ui + vj 

Equation (3) may be expanded, after substitution, 
to yield; 

“ {[z(v+Q^u) + v(n^v-JlyU)-z^ ft^]i 


Since small deflections are assumed for the 
beam, the reflector can be assumed to be located 
at a constant distance from G, the Shuttle mass 
center. 

Using the transfer theorem fjgr the angular 
momentum^, the angular momentum, of the 

reflector, assumed rigid, about G can be expressed 

\/G “ ^R/Gj \/Rq “r * dt I Rq 

where Ik/Gi* ^ke inertia tensor of the reflector 
^xpres^ed it G^, its center of mass, and Gr/Pq “ 
^R/S^ reflector’s initial 

angular velocity, its angular velocity relative. to 
the Shuttle, and the Orbiter’s inertial angular ve- 
locity) are both expressed in the same coordinate 
system, R 2 , moving with the reflector. 

D. Angular Momentum of the System About G 


+ [- z(u -G v) + u(a u-fi v) + z^« ]j 
z y X y 

+ [u(v+ Ji^u) - v(u-Q^v) *f z(uG^+vGy)]k} dm (4) 
where 

u(z,t) • I p^(t) (z) and v(z,t) « 

n ^ 

I p^(t) S^(z) 
n ^ 


The total angular momentum of the mast about 
G is obtained by integrating Equation (4) over 
the total length of the mast. 


“h/o- /"■ ■*./<: 

To simplify the notations, let 


r- /e%\ — f A /▼ cosBL slnBLix 
f^(e) - iA^a -g 


„ sinBL . cosBL 1^ 
+ B, (L - Jt) 


M r f ^^kBL 


» coshBL X 

L, - ) 


+ D, (L 


sinhBL 

B 


coshBL ,1 m 
__ + ^)} 


After substitution of f^ and — for p, ^here 
M * mass of the mast, in the expression of 
and considering only the effect of a single modd 
(for demonstration purposes here, with frequency 
«) one arrives at; 


^ The angular momentum of the system about G, 
Hsyst/G» given by the sum of the angular momen- 
tum or each of the three components evaluated about 
the same point, G. 

III. Rotational Equations of Motion (Torque Free) 

The rotational equations of motion for the 
system, when free of all external torques, are 
obtained as; 

* • 

^^syst/G^I^O “ ®syst/G/S * ®syst/G “ ^ 

(9) 

The vector equation (9) Itself is equivalent to 


H+GH-OH- 
x y z z y 

0 

(roll) 


H+GH-OH- 
y z x X z 

0 

(pitch) 


+ n H - G H - 
z X y y X 

0 

(yaw) 

(10) 


When the gravity-gradient torques taken at the Shut- 
tle center of mass are included in the equations 
of motion, the linear system dynamics appear in the 
following state former 

X « AX -I- C (11) 

where C is a constant and X* (V' ,9 ,6 , involv- 

ing the Euler roll, pitch, and yaw angles (and their 
rates) , respectively. 
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This indicates that the system equilibrium posi- 
Mon is no longer: ’Pn " ®0 “ *0 

offset (X,Y) in the attachment point of the beam 
to the reflector (Fig. 1). 

Let , 6 , and ^ be the equilibrium posi- 
tion for tSis configuration of the system. Then, 


^ • 

♦e+ 

and « 

n 


e * 


; e - 

h2 


^ • 


; 1 - 

^3 


The new 

state vector is [n^ 

? ^2 ’ ^3 ’ 


Also 

6 , and ^ 
e ^e 

satisfy 



*1 

K + ^2 

e + 

e 

^3 ' 

■®19 

*7 


0 + 
e 

*9 *e ■ 

■*20 

®13 

K * *14 

6 + 
e 

*15 *e * 

-*21 


where is the mode shape vector; L is a 
1 liif.'ir opor/itor, which wli£n applied to q yields 
the elastic force on dm, f represents the gravita- 
tional force per unit mass, and *e represents the 
external and control forces on dm. 


After substitution of the values for the in- 
tegrals into Equation (12) and rearrangement of 
the terms, the generic equations are obtained 

in the following form:^’® 


A + 0) A + ♦ /M + I ♦ /M 
n nn nn . mnn 
m«l 


= [g„ + £ g + E ]/M 

®mn n"* n 

m=l 


(13) 


where A^, ui , and M are the n^^ modal anq>li- 
tude, frequency, shape functions and modal 
mass, respectively. 

the effect of gravity on the n^h mode 

the effect of control and external forces 
on the n^^ node 


this simulataneous system is solved using 
[a] - [a] e^, ®e’ “ 

Ca"^] ta] 


®ran 


the coupling effect of gravity from the 
mode on the n^^ mode 


_ the coupling effect of elastic forces 
ton 

from the m^^ mode on the n^^ mode 


After substituting the new state vector in 
the equations describing the system dynamics, 
linearizing them about the new equilibrium posi- 
tion, recasting them into state format, one 
arrives at a system which can be cast in the 


- 

A’n 

(where a^ 

- *is 

are constants)^ 

0 

0 

0 

1 

0 

0 


0 

0 

0 

0 

1 

0 


0 

0 

0 

0 

0 

1 


*i 

*2 

*3 

*4 

*5 

*6 



*8 

*9 

*io 

*11 

*12 


*13 

*i4 

*l5 

*ie 

*l7 

*18 



The open-loop system in this configuration is 
unstable due to the unfavorable Inertia distribu- 
tion. 

IV. Generic Mode Equations 

The generic mode equations are obtained by 
taking the modal components of all internal, exter- 
nal and Inertial forces acting on the system, i.e., 

•• • * 

/ + r + 2o)xr + o)xr + a)x(wxr) dm • 

M _ _ 

/ • CUq)/^ + f +i] dm (12) 

M 


V. Control of the Orbiting SCOLE with . 
the First Four Modes Included 

This model of the SCOLE is assumed to be con- 
trolled through the three torquers on the Shuttle 
and the six actuators located by pairs at z*-L/3; 
z *~2L/3 on the mast: and at G. , the mass center of 
the reflector (Fig, 1). The pairs of actuators are 
arranged in such a manner that one acts along the x 
direction and the other in the y direction. The 
actuators, when activated to provide vibration con- 
trol to the mast, will develop torques about the 
Orbiter center of mass. Each actuator provides a 
maximum of “ 800 Ib.^ force; the resulting 

torque contr^utea by all six actuators is computed 
as 

= V 2 V2y/3 + V3y) 1 

- <"1x^3 + 2 v2^/3 + - (YF, 

- XF V3y)k (14) 

This is added to the torques provided by the^ Shut- 
tle’s three torquers: T^-M G,.i + H^Uy3 

+ M U k, where M ^ 10, 000ft. -lb. ^ to 

yield the total available control torque for the 
system as: 

t » [M^ ^ <^ly/3 + 2v2y/3 

+ v^y)]L[My Uy - F^L(v,^/3 + 2v2^/3 
+ V3^)]j‘ + [M^ U, + X Fy V3y - Yv 3^ F^]k (15) 
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AIAA COPY SHEETS 


with the control and global state vectors, res- 
pectively, chosen as 


U 




'2x’ 


'2y* 


V 


3x* 


3y 


II , II , II 

X y /. 


r • # * 

and X " > ^2* ^3* ^1* ^2* ^3* ^4* ^1* ^2* ^3* 

• • • • T T 

A » A2 f , A^ J 

with l|vjjjll< 1; llvj^ ||< 1; and ||U Jill, v 
represents the force due to the ith actuator in 
the X direction and similarly for v , i « 1,2,3. 
The control influence and the system^ state matri- 
ces are obtained, respectively as: 


B* 

14x9 


A 

14x14 


0 

» 

0 


7x6 

1 



B^(F,L,X,Y) 

1 ®2(«x*»y>“z) 

3x3 

®3 4x6 

1 

4x3 


0 i 0 1 

I. 1 o’ 



--1- — r 
0 1 0 

— 1 

°_L !* 







^5 ! ^ 1 

A- 1 0 

‘ 1 




(16) 


(17) 


The specific submatrices are explicitly defined in 
terms of the SCOLE parameters in Ref. 5. 


The largest disturbance in the flexible nodes,, 
caused by an attitude variation is observed dur- 
ing the roll axis maneuver (Figs. 2-3); the first 
mode 1s the most excited; its amplitude doesn't 
exceed 0.13 ft. (0.12 of L) . All the transients 
are damped out within 23 seconds, largely due to 
the contribution of the additional 2 pairs of 
actuators located on the mast at * -L/3 and z ■ 
-2L/3, During this response the reflector **y'* axis 
actuator provides a maximum of 21C lb. while the 
forces in the two "y” actuators located a z ■ -2L/3 
aid z*-L/3, reach 120 to 80 lb, respectively. The 
Shuttle, *'x" torquer provides a maximum of 2800ft. -* 
lb torque bringing to 52,500 ft. -lb the maximum 
value of the x component of the composite control 
torque required for this maneuver* 

In turn, the first flexible mode was given 
an initial amplitude equal to 1.0% of L, to stay 
within the linear range, Figsi.6-7 show the tran- 
sient responses. Also depicted is the result of 
intraflexible modal coupling (Fig. 6) • For this 
control strategy, the disturbance in the flexible 
mode, for the initial conditions considered herein, 
is damped in 15 seconds while its effect on the * 
attitude angles takes almost 25 seconds to dis- 
appear . 


VI * Conclusions 


1. The SCOLE system, with gravity-gradient tor- 
ques included in its open-loop dynamics, is 
unstable. This is due to the inertia distri- 
bution of the system in the particular con- 
figuration where the Shuttle roll axis ncminalljr 
follows the orbit. 


Here, a control, U, which minimizes the per- 
formance index 

J - / (x’^QX + U^RU) dt is obtained after 

7 

using the ORACLS' package to solve the steady statei 
Riccati matrix equation* 

The equations describing the closed-loop sys- 
tems, X - A'X + BU have been numerically inte- 
grated and the corresponding mathematical model 
simulated for 

Q - diag. [5x10^,5x10^,5x10^,5x10^,5x10*, 

5x10^ ,5x10^,10,10,10,10,10,10,10] 
and R as diag. [10,10,10,10,10,10,1,1,1] 

Since both position and rate feedback of the Shut- 
tle rotational motion and beam elastic motion will 
be utilized it is logical to place a greater 
penalty on the position displacements. Also, 
since the roll (and to some extent the pitch) are 
easier to excite than some of the elastic motions 
it seems intuitively correct to relax the penalty 
of the Shuttle control inputs as contrasted with 
the remaining control penalty elements. 

The transient responses to some Initial per- 
turbations, confirm the controllability of the 
flexible SCOLE system. During the simulation of 
this model, the three variational attitude angles 
(roll, pitch, and yaw) are each subjected to a 6 ^ 
single axis displacement. For each case, the 
effects of such displacements on the modal ampli- 
tudes of the first four modes are studied. 


'2. In the absence of control forces and torques, 
the system will oscillate about an equilibrium 
position biased from the Shuttle's symmetry 
axes. This phenomenon is due principally to 
the offset of the beam attachment to the re- 
flector from the reflector's mass center. 

3. A robust control law based on the linear regu- 
lator theory can be implemented for station 
keeping with maximum control efforts signifi- 
cantly below saturation levels. 
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Fig. 7 Linear Model of SCOLE with Flexibility 
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RAPID SLEWING OF THE ORBITING 
SPACECRAFT CONTROL LABORATORY EXPERIMENT 
(SCOLE) USING LQR TECHNIQUES 


Abstract 

The rotational equations of motion, describing 
the dynamics of the (rigidized) proposed orbiting 
Spacecraft Control Laboratory Experiment during 
the station keeping phase, are derived using the 
Eulerian formulation. When the attitude angles 
(roll, pitch, and yaw) are assumed small, a sta- 
bility analysis is conducted for the system. It is 
seen that the pitch equation decouples from the 
roll and yaw equations when the interface between the 
mast on the reflector is not offset or the offset 
is only along the Shuttle roll axis. When a 
second offset is introduced along the pitch axis 
the system and when the gravity-gradient torques 
are present in the dynamics, the system assumes a 
new equilibrium position. The linear regulator 
theory is used to derive a control law for the li- 
tnear model of the rigidized SCOLE. This law is 
applied to the nonlinear model of the same con- 
figuration of the system and preliminary single 
axis slewing maneuvers (20o amplitude) are simu- 
lated. 

I. Introduction 

In several versions of a recent paper, 
jSCOLE(l) (Spacecraft Control Laboratory Experiment), 
[Lawrence W. Taylor, Jr. and A.V. Balakrishnan have 
)d escribed the ground based experimental program 
Iwhich is to be used to test techniques derived for 
active control of flexible space structures. It 
is a laboratory experiment based on a model of the 
Shuttle connected to a flexible beam with a re- 
flecting grillage mounted at the end of the beam 
(Figure 1). The interface connecting point between 
[the reflector grillage and the beam is offset in 
jtwo cartesian directions with respect to the center 
pf mass of the reflector. As a part of the design 
challenge, the authors stressed the need to di- 
jrectly compare competing control design techniques 
pnd discussed the feasibility of the such direct 
comparison. The challenge consists in slewing the 
laforementloned model line of sight through 20^ in 
tiinlmum time. 

Based on the equations describing the motion 
pf the SCOLE system, provided in reference 1, the 
[expression for the reflector line of sight (LOS) 
terror was expanded analytically and studied 
parefully. The results showed that the SCOLE *s 


LOS error is independent of the Euler yaw attitude 
angle so, only two, instead of originally three, 
angular parameters were needed to be concerned 
with in designing the pointing .slew, maneuvexe. It 
was also suggested that a two stage control strat- 
egy which would first slew the whole system, as if 
rigid and then damp out the residual undesired 
mast vibrations, is most appealing. The numerical 
simulation test results of Reference 2 indicated 
that the single axis bang-bang or bang-pause-bang 
slew maneuvers work fairly well for pointing the 
LOS of SCOLE. The best pointing accuracy and 
shortest slew time were attained when using the 
Shuttle torques and the actuators placed on the 
reflector while imposing a 5 degree/second slew 
rate limit on the design. 

In the present study, a mathematical model of 
the SCOLE system is developed assuming: the 

Shuttle, the mast, and the reflector to all be 
rigid. This development is based on the Eulerian 
approach. The technique consists in isolating an 
elemental mass of the system and deriving Its 
angular momenttan taken at the mass center of the 
0rbiter.(3) The expressions obtained for the 
elemental masses of the components of the system 
are Integrated over the mass of the entire system 
to yield its angular momentum about the center of 
the Orbiter. The time derivative of the system 
angular momentum is then equated to the gravity- 
gradient torques^^) of the system about the same 
point. Such a vectorial equation, when projected 
along the three axes of rotation, yields the sys- 
tem rotational equations of motion. These equa- 
tions are then linearized to yield a model which 
provides the basis for the control law synthesis 
developed in this study. 

II. Development of Equations of Motion 

The rotational equations of motion for the 
system will be derived by taking the time deriva- 
tive of the angular momentum of the system at G, 
the center of mass of the Shuttle, and by equating 
it to the external torques applied to the system. 

A. Angular Momenttsn of the Shuttle about its 
Mass Center, G. 

The angular momentum of the Shuttle, taken as 
a rigid body, about its center of mass, G is 

®S/G “ ”s/G S/Rq 

where I . is the inertly tensor, with elements 
Igi, of the Shuttle and the Shuttle's 

Inertial angular velocity, with components 
i * x,y,z, along the Shuttle's symmetry axes. 
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B, Anuglar Momentum of the Rigid Beam about G 

The angular momentum of the beam about the 
Shuttle mass center » G» can be expressed as 

■ Vg " “m/g \ “m/g 

tensor of the beam transferred at G, using the 
parallel axis theorem. 

Ml 2 


«M/G “ - T ^ ^ S 


( 2 ) 


vhere M is the mass of the beam whose length is L. 
C. Angular Momentum of the Reflector about G 

Since here the beam is considered rigid, the 
reflector can be assumed to be located at a con- 
stant distance from G, 

3 

Using the transfer theorem ^for the angular 
momentum, the angular momentum, /«, of the 
reflector, assumed rigid, about G can be ex- 
pressed as: 


H, 


d 


R/G 




(3) 


where » the* inertia tensor of the reflector 

jexpressed ia£ Gj^, its center of mass, and « 

^S/R inertial angular ^ 

velocity of the reflector. 

Angular Momentum of the System about G 

^ The angular momentum of the system about 
G,Hgyst/G» is given by the sum of the angular 
momentum of each of the three components evaluated: 
about the same point, G: 


^syst/G *' ^S/G *Sl/G \/G 
®syst/G ■ ^ V ^ 


(4) 


E. Rotational Equations of Motion (Torque Free) 

The rotational equations of motion for the 
system, when free of all external torques, are 
obtained as: 


dt ^”syst/G'|RQ » ”syst/G/R 
^S/Rq * ®syst/G " ^ 


(5) 


The vector equation (5) itself is equivalent to 


x y z z y 


H+HH-GH-O 
y z X X z 


H+nH-«H*0 
z X y y X 


( 6 ) 


III. Stability of the SCOLE 
in some of its Configurations 

Equations (6) describe the torque-free non- 
linear dynamics of the SCOLE configuration. In 
what follows, the stability analysis of the 
"rigid ized" SCOLE system will be conducted in threes 
different steps for the cases where the Shuttle's 
roll, pitch, and yaw displacement amplitudes 
respectively) are assumed small. 

First, it will be assumed that the interface 
point between the beam and the reflector is the 
reflector center of mass; second, still assuming^ 
IChe mast‘rlgid, the interface point will be offset 
in the "X" direction; finally, a two dimensional 
offset of the interface point will be introduced. 
The mast will still be assumed rigid. The system 
dynamics, in all the aforementioned cases, includes 
the gravity-gradient torques.^ 

A. The SCOLE System without Offset . 

In the absence of offset in the location of 
the interface point (X**Y»0) , with gravity-gradient 
torques present in the system. Equations (6) can 
be rewritten as^: 




- I, 


+ - I, 


ML^/3 + 3(13 - Ij)] - 0 ( 7 ) 


e [lg^+Ijj^+Mj^LW^/3]+3wJe(I^-l3)+3oQl^- 0 (8) 

4 I + ♦ (I- + 1r ) + 

^ ^83- ^83 - ^R 3 - ^R 3 - ^R3^ - h . 


+ - I. 1 + 3lJ 


(9) 


where is the Shuttle's xz product of inertia 
and 2 3 composite system inertias with 
respect to the Shuttle axes. It is seen that in 
, such a configuration, in the linear range, the 
: equation describing the pitch motion (Equation 
(8)) of the system decouples from the equations 
describing the motion in the two remaining degrees 
of freedom (Equations (7) and (9)). 


Equation (8) can be recast in the following 


form : 


6h^ - 0h2 + h^ 


( 10 ) 


in which, 
h 


>s, ^ s, - * ¥ 


h2 ' 3u)q (I^ - 1^) and h^ 
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The homogeneous part of Equation (10) yields the + - I. + I- + I_ + - I ] 

following solution: ^ ^1 ^2 ^3 *1 


Cie 


6t 


C2e 


-6t 


where 


6 = /h^7h^ 


since for this configuration, h^/h^ > 0, Qj^(t) 

is unstable. Thus, in the absence'^of control, 
the system is seen to be unstable in its pitch 
degree of freedom. 


Equations (7) and (9) which have the follow- 
ing forms, respectively. 


- “o« v'l 

- Uq i|. {3(Ig + Mj^XL)} - 0 (13) 

4 

where I. represent the composite inertias about 
the Shuttle axes for this case and Wq is the 
(circular) orbital angular velocity. 


Again, it is seen that in this configuration, 
the pitch equation, (Equation 12 >, decouples from 
the roll, (Equation II)) and yaw (Equation 13)) 
equations euid can be rewritten as: 


+ ^V.2 • - ’l^k^ * 0 

.. .. ^ 

+ ^n^ - + \|;n^ ■ 0 

can be recast in the following state matrix 
format ; 




0 0. 1 0 





0001 





P2 P4 Pi P3 


• 



.-P 6 -Ps -P5 -P) . 


1 

•0^ 


Some of the eigenvalues of the state matrix, 
in this subcase have positive real parts, based 
on the actual SCOLE system parameters indicating 
instability in the open loop dynamics of the roll 
and yaw degrees of freedom. 


eh| - eh^ + - 0 (14) 

where, - Ig + Ij^ + Mg(X^+ L^)+ ML^/3 

‘‘a ■ 

4 

+ Mj^XL) 

Here again, b2/h. is a positive quantity. By " 
analogy with the previous configuration, 

a 

9(t) “ cosh 6’t sinh 6*t *f h-/h^ 

U gf j ^ 

(1-cosh 6*t) (15) 


with 


6 ’ 




In the absence of control, it is seen that the 
pitch angle is unbounded indicating an instability 
in that degree of freedom. 


B. The SCOLE System with Offset in the *'X '* 
Direction 


The configuration analyzed in the previous 
section .is upgraded to the one considered here 
by letting X be non-zero in the equations of 
motion (6) and by setting the "Y offset" equal 
to zero. 


The equations of motion then become: 

Vtlg + 

-'♦<'54 + "e”'> - ♦•ot'Sj * ^3 - 's, 

*\ * \ - V"’ ■ «S*''S3 

- 'S, - ^ -V* *^"3 - ■ ” 

( 11 ) 

+ 3u)J(lg + M^XL) - 0 (12) 

4 

4 Clg + MjjXL] +'i [Ig + Ir + 


A reasoning similar to the one previously 
done for the case without offset, enables one 
to recast Equations (11) and (13) in the following 
state matrix format: 



(16) 


I ^ t 

p., 1 « 1 8 are defined in terms of the k^ andn^ 

similar to the case without offset, and appro- 
priate kj^, n^ now Include the effects of the X 
offset.^^ Here again it is seen that some of the 
eigenvalues of the state matrix. A*, have positive 
real parts. Therefore, the open loop dynamics of 
the system are seen to be unstable In its roll 
and yaw degrees of freedom. 


C. The SCOLE System with Offset in Both the "X" 
and **Y" Directions 


If once more the description of the system 
dynamics is upgraded by introducing the **Y offset” 
the rotational equations of motion become: 


3.3 


OF POOR QOALiTY 


T '''>1 -♦'xz 

-eM^XY -;u.oCis^+ IS3 Ir^ Ir^ 

-“o®V" - “ 0'^t‘s3- IS3^ 'R3 - Ir. 

- r ^ 3('zz- ^y>3-“0*^xz-^®' 

4u,2[M^yL + 31 ] - 0 


xy 

(17) 


e[l5^+ + \ + Mj^a^ + X^)] + 

•h|<MjjXy + o)q*Mj^XY - ojq^MjYL + Ugf^YL - 3 u.qU.I 

+30)^6 (I_-I..) + 3o)„I 


XX zz 


0x2 


xy 

(18) 


•Hsj* h* "r“"* '■'>1 - «zx 

%- \* \* \- \ 

+2M^Y^] - 9MjjYL - 0)^6 Mj^XY + 30)^61 

-“0*«-»> - 4*\- V '«i‘ "'>: 

-H1r(V^-X^)3 - »^\XY - 0 


yz 

■2 


(19) 


It should be noted here that the pitch equa- 
tion no longer decouples from the roll and yaw 
equations. Equations (17), (18), (19) can be 
recast in the following state matrix format 

X « A** X + C or 


. .‘^*'1 new state vector Is [n,, n,, n,, n, , 

n2> ^3]^* Also e^, and satisfy ^ ^ 


®1 ♦e ®2 ®e + *3 ♦ 


e “l9 


'*7 "^e •" “8 ®e ®9 *e ' ®20 


^13^ \4®e ■*■ \5*e ' *21 

and this simultaneous system can be solved to de- 
termine 0g, 


After substituting the new state vector in 
the equations describing xhe system dynamics » 
linearizing them about the new equilibrium posi- 
tion, recasting them into a state format , one 
arrives at [n.j * tA ^ [n^] . The real parts 
of three of the new state matrix eigenvalues are 
found here to be positive indicating that the open 
loop system in this configuration is also unstable* 


IV. Control Synthesis 

First, within the linear range, the motion of 
the rigldized SCOLE is controlled using a strategy^ 
based on the linear regulator problem when the 
system is subjected to some small perturbations in 
Its degrees of freedom; second, the control stra- 
tegy derived for the linear model of the rigldized 
SCOLE is applied to the non-linear model of the 
same configuration. Preliminary slew maneuvers 
are tested by assuming single axis initial per- 
turbations of 20® in the roll, pitch, and yaw de- 
grees of freedom, respectively. The three Shuttle 
torquers and the two actuators on the reflector 
(Fig, 1) are then assumed to be the only sources 
of control moments. The controllers are seen not 
to reach saturation. 
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.^2lJ 


( 20 ) 

where, the a. are functions of the various com- 
ponent inertias and the X,Y offset parameters. 5 


Since the Shuttle axes do not correspond to 
the principal axes of the system, the system 
dynamics appear in the following state form: 


X “ A .. ^ 

offset 


X + C 


( 21 ) 


where A^^^^ ■ A*' form Eq. (20) indicating that 

the system equilibrium position is no longer 

'^'0 “ ®0 " *0 ■ 


A. Control of the Linearized Model of SCOLE 


During the control of this model, it is 
assumed that the actuators located on the mast 
(proof masses) are not activated. As a result, thet 
system is controlled by means of the Orbiter 
torquers and the actuators located on the re- 
flector (Fig. 1), 

Since the Shuttle is equipped with three 
torquers acting about the x,y, and z directions, 
the total control torque available can be written 
as 

T + 130Fy Vy) i + (MyUy-130 v^)j 

+ (M^U^ + 32.5F^v^ + 18-75Fy Vy)k) ft. -lb 

with the limits for M , M and M - 10,000. ft. 
lb; F and F • 800 Ib.^ The constraints, there- 
fore, *are ^ 

lUylii; lUjIli; 

Ivyli 1 


Let ij;^, e , and 4> be the equilibrium position where U, the control vector is expressed as 
for this configuration of the system. Then, « 

U « [vj^, v^, U-, U^, while the control 

\i> • ^ influence mat rix^can^ then be written as: 

e • + ^2 0 * 02 

^3 ^ * ^3 
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0 

0 

0 

0 

-130F^ 
32. 5F.. 


0 

0 

0 

130F 

y 

0 

18.75F 


The optimal control U which minimizes a per- 
formance index 


J - f (x'^QX+u'^RU) dt 

0 

is given by 

U - -KX - -(r"Vp)X 

where P is the positive definite solution of the 
steady state Ricatti matrix equation.^ The 
equations describing the closed loop system can 
be recast in the following matrix format: 

X - AX + BU 

where A = and X = h 2 > nj]'^ 

from the discussion following Eq. (20). After 
substitution of -KX for U, the closed loop equa- 
tion can be rewritten as 

X • (A-BK)X 

A parametric study was conducted by first 
examining the variation of the real part of the 
least damped mode as a function of different 
values for the (assumed) diagonal Q and R 
weighting elements (Figure 2), In this initial 
study, each of the diagonal Q elements were 
assumed equal i.e. Q»diag. [SQ] and also each of 
the diagonal R elements were assumed equal R « 
diag. [SR]. Figure 2 corresponds to a model of 
the rigidized SCOLE system where the dimension- 
ality of the state vector is 6 x 1 and 3 Shuttle 
torquers plus 2 reflector actuators describe the 
control Inputs. 

It can be seen from Fig. 2 that the best 
closed-loop transient results are obtained from 
using larger values of the state penalty along 
with smaller values of the control penalty ele- 
ments. However, when the closed loop dynamic 
responses were simulated using the best combina- 
tions of Q and R it was seen that some of the 
controllers reached saturation levels for re- 
sponses with initial conditions on pitch, roll, 
and yaw taken within the slewing angle range 
(i.e. approx. 0.3 rad.). 

As an alternative, the concept of split 
weighting of both the state and control penalty 
elements was considered, initially for the rigid- 
ized SCOLE model. Since the roll (and to some 
extent also the pitch) are easier motions to 
excite than the yaw, due to the SCOLE moment of 
Inertia distribution, it seems intuitively correct 
to relax the penalty of these control Inputs as 
contrasted with the remaining control penalty 
elements. Also since both position and rate 
feedback of the Shuttle rotational motion will be 
utilized, it appears logical to place a far greater 


penalty on the (angular) position displacements. 
Based on this philosophy and by trial and error, 
the set of Q and R which produced the largest 
absolute value of the real part of the least 
damped mode (while at the same time avoiding sat- 
uration during 20 <^ single axis slewing maneuvers) 
was selected as 

Q = diag. [5x10^2, 5x10^^, 5x10^^, 1, 1, 1 ] 
and R * diag. [l, 1, .1, .2, l] 


For this set of Q and R the closed log eigenvalues 
for the rigidized SCOLE model ar^ Q^lCjulated 


R(xp 

-0.431436E+02 

-0.431436E+02 

-0.132023E+03 

-0.132023E+03 

-0.328320E+03 

-0.328320E+03 


Im(X^) 

0.431436E+02 

-0.431436E+02 

0.132023E+03 

-0.132023E+03 

0.328320E+03 

-0.328320E+03 


It has been assumed here that all the state vari- 
ables are available at each instant (observabi- 
lity matrix ■ ^ 5 )* 

The closed loop dynamics has been simulated 
for transient responses to a 6 ° initial perturba- 
tion in roll. Figure 3 shows that a 6 ° perturaba- 
tion in roll is damped out in approximately 13 
seconds. During that single axis maneuver, it 
should also be noticed that the coupling disturbs 
the yaw degree of freedom, which reaches a max- 
imum amplitude of 0.25® degree. Figures 3a, 3b, 
and 3c show, for the 6 ® maneuver about the roll 
axis, the forces required from the reflector 
actuators, the efforts produced by the Shuttle's 
torquers, and the components of the equivalent 
total torque acting on the SCOLE system, respec- 
tively. The reflector *'y** actuator and the 
Shuttle's "xV torquer are the more active con- 
trollers for this maneuver, as expected. 


B. Rigidized SCOLE Preliminary Slew Maneuvers 

In this section, the equations governing 
the motion of the rigidized SCOLE, outside of the 
linear range, are developed from the most general 
rotational equations of motion previously derived. 
The control laws obtained from the application of 
the linear regulator theory to the linearized mo- 
del of the rigidized SCOLE are tested for large 
amplitude manuevers. The closed loop system 
dynamics are numerically simulated. For single 
axis slew maneuvers about the roll, pitch, and 
yaw axes, respectively, the time responses for 
the Euler angles, the control efforts required 
of the reflector actuators, control torques 
demanded from the Shuttle's torquers, and the 
components of the total control moments, are 
depicted in the subsequent figures* This enables 
one to determine the margin left in which to opti- 
mize the control strategy without causing satura- 
tion of the controllers. 


The equations governing the motion of the 
rigidized SCOLE system during large amplitude 
maneuvers in the presence of gravity-gradient 
and control torques are obtained as: 


OF POOR QUALITY 
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±) The Roll Equation 


«x + «y «z - »y “ 


ii) The Pitch Equation 
H^+ n H - n H - 

y z X X z 

iii) The Yaw Equation 


z X y y X z 


( 22 ) 


(23) 

( 24 ) 


where T , T , and T^ are the components of the 
external torques acting on the system (including 
the control torques previously derived for the 
linear model of the rigidized SCOLE where the 
feedback, now depends on the original Euler angles 
and their rates for maneuvers made relative to 
the Shuttle roll, pitch, and yaw axes) , 


The closed-loop system dynamics described by 
Equations (22), (23), and (24) have been numer- 
ically simulated and the results are shown in 
Figs. 4 to 7. Fig 4 shows the time responses to an 
initial 2(P alignment in roll. It is seen that a 
20® slew about the roll axis can be achieved in 
about 30 seconds. 


For this control strategy, the single axis 
slew maneuver about the roll axis uses 80% of the 
control forces available from the corresponding 
actuator located on the reflector, and 80% of the 
control torque available from the corresponding 
Shuttle torquer. None of the controllers reach 
saturation. Also depicted in the Fig. 7 are 
the components of the total control moments for 
this case (moments of the reflector control forces 
taken about the Shuttle’s mass center, plus moments 
of the Shuttle’s torquers) . This will make pos- 
sible a comparison between this strategy and other 
future control laws which would be based on the 
two point boundary-value problem,^ where this 
or combinations of control inputs may be employed. 
Similar results have also' been simulated for slew 
maneuvers about the Shuttle’s pitch and yaw axes.^ 

V. Conclusions 


In conclusion, it is seen that: a) a control 

strategy derived from the linearized model of the 
rigidized SCOLE based on the linear regulator 
theory, works well when It Is used for single 
axis slew maneuvers through amplitude angles 
as large as 20®. For the case of the yaw axis 
slew maneiiver^ (wh^re the control effort does not 
benefit directly from the long moment arm provided 
by the mast) by suitable adjustment of the gains. 

It Is possible to achieve a slew in approximately 
100 secs, with the use of up to 90% of the Shut- 
tle’s (yaw) saturation torque. In comparison 
with the results of Ref. 7, where fast slews of 
order of 10-20 secs are accompanied by series of 
bang-bang type control efforts In more than one 
control input, the results here provide a definite 
trade-off between the slightly increased slew times; 
‘with the considerable reduction in over-all con- 
trol effort; b) the SCOLE system, with gravity- 
gradient torques included in its open-loop dynamics, 
is unstable. This is due to the inertia distri- 
bution in the particular configuration where the 
Shuttle roll axis nominally follows the orbit; 
and c) the equation describing the pitch motion 


decouples, within the linear range, from the roll 
and yaw equations, when the gravity-gradient tor- 
ques effects are present In the system dynamics, 
and when the mast attachment point on the reflec- 
tor is not offset or when the offset is parallel 
to the roll axis. 
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FIGURE 1: DRAWING OF THE SCOLE CONFIGURATION 


FIGURE 2: REAL PART OF LEAST DAMPED MODE VS. [R] 

ABS.(R X) Absolute Value of Real Part 
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A NUMERICAL APPROACH FOR SOLVING RIGID SPACECRAFT 

MINIMUM TIME ATTITUDE MANEUVERS 
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Abstract 

The minimum time attitude slewing motion of a rigid 
spacecraft with its controls provided by bounded torques and 
forces is considered. Instead of the slewing time, an integral 
of a quadratic function of the controls is used as the cost 
function. This enables us to deal with the singular and 
nonsingular problems in a unified way. The minimum time is 
determined by sequentially shortening the slewing time. The 
two-point boundary-value problem is derived by applying the 
Pontryagin's Maximum Principle to the system and solved by using 
a quasilinearization algorithm. A set of methods based on the 
Euler's principal axis rotation are developed to estimate the 
unknown initial costates and the minimum slewing time as well as 
to generate the nominal solutions for starting this algorithm. 
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It is shown that one of the four initial costates associated 
with the quaternions can be arbitrarily selected without 
affecting the optimal controls and^ thuS/ simplifying the 
computation. Several numerical tests are presented to show the 
applications of these methods. 

Introduction 

The problems of large-angle attitude maneuvers of a 
spacecraft have gained much considerations in recent years. 

In these researches / the configurations of the spacecraft 
considered are: (1) completely rigid/ (2) a combination of rigid 
and flexible parts, or (3) gyros tat-type systems. The 
performance indices usually include minimum torque integration, 
power criterion, and frequency-shaped cost functionals, etc. 

Also some of these investigations utilize feedback control 
techniques. In this paper, the minimum time attitude slewing 
control problem of a rigid spacecraft is considered. 

In Ref. 2, the problem of the rapid torque-limited slewing 
of the rigidized SCOLE^ about a single axis (x-axis) is 
considered. The control torque about this axis is of a bang-bang 
type or a bang-pause-bang type. The control laws are developed 
based on a simplified model of the SCOLE and then used on the 
practical model (with nonzero products of inertia); hence, this 
leads to a large error of the attitude after the slewing. Also 
it seems that no details were given for the controls about the 
other two axes (y, z). 
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In the present paper, the optimal control theory (Maximum 
Principle) is applied to the slewing motion of a general rigid 
spacecraft (including the rigidized SCOLE, without 
simplification). The slewing motion need not be restricted to a 
single-axis slewing. The computational procedure based on a 
quasilinearization algorithm is developed to solve the resulting 
two-point boundary-value problem. 


Euler Rotation and State Equations 


The attitude of a rigid spacecraft can be described by 
either a quaternion vector q=[qo <31 <32 < 33 1 ^/ which satisfies a 
constraint equation, q^q=l, or a direction cosine matrix C, 

ql+ql-ql-ql 2 (qiq2+qo‘53 ) 2 (qiq3-qo<32 ) 

2(qiq2-<5o‘^3) 2(q2‘53+‘^o‘^l) 

2(qiq3+qo‘32) 2(q2‘i3“‘3o'3i) 

It is known that a single axis rotation of the spacecraft 
can also be expressed by a quaternion 

qo=cos(0/2); qj[=eiSin(0/2 ) , i=l,2,3 (2) 

where 0 is the angle of rotation , and e^ are the direction 
cosines of the rotation axis. 

The Euler rotation theorem tells us that an arbitrary 
orientation of a rigid body can be realized by rotating it about 
a principal axis (eigenaxis) through a certain angle from its 
initial position. The desired rotation quaternion, q, between 
the initial position q(0) and the final orientation q(tf) can be 
obtained by the following equation 
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qo ■ 


qoo qio q 2 o qio 


qof 

qi 


-qio qoo q3o -^20 


qif 

q2 


-q 2 o -q3o qoo qio 


q2f 

. ^3 . 


-q3o q 2 o ~qio qoo . 


. ^3f . 


where the second subscript "q" "f " represent the Initial 

time and final time, respectively. The angle of rotation and the 
unit vector, e=[ex ©2 ©sl^r along this eigenaxis are 

0*=2cos"iqo; ei*<3i//l-qo^» i=l,2,3 (4) 

The equations of motion of a rigid spacecraft are 


4=(l/2)Sq (5) 

I u =c3icj+bu ( 6 ) 

where 6 J=[Wx CJ 2 ^ 3 ]^ is the angular velocity vector, B is a 
3xn control influence matrix, u=[ui U 2 U 3 ••• u^l^ is the 
control torque and force vector, and 
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Premultiplied by the inverse of I, Eq. ( 6 ) can be rewritten as 


cj =I-1 Ci6J+i-1bu (7) 

The boundary conditions for the states, q and 6 i, are given as 

q( 0 ), CJ( 0 ); q(tf),CJ(tf) ( 8 ) 
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optimal Control 

Conventionally f the time optimal problem involved here is 
to seek a solution of Eqs. (5,7) satisfying the boundary 
conditions (8) and minimizing the slewing time 

tf» r (l)dt 
0 

The control variables must satisfy the constraints 

luj^uimax“'^ib» i~l»2, •••,n. (9) 

The Hamiltonion for this problem can be written as 
H=l+(l/2)pTC q+rT(i-lClC*)+i-lBu) 
where p and r are the costate vectors associated with q and U, 
respectively. They satisfy the necessary conditions: 

p=-(3H/3q), r=-OH/3«0) (10) 

By Pontryagin's Maximum Principle, the optimal control 
minimizing H can be determined by 

Ui=-Uij 3 sgn(B^I“^r )i, i«l,2, ..., n. (11) 

which are of the bang-bang type, except for the singular case in 
which (B^I“lr)i=0 over some .non- zero time intervals. For the 
nonsingular bang-bang optimal control problem, a shooting 
method^O was tried and failed due to the nonexistence of the 
inverse of a partial differential matrix. 

The singular problem does occur if a special case is 
considered in the following. When B is a 3x3 unit matrix and I 
is a diagonal inertia matrix, the control in Eq. (11) is 
simplified as 

u^=— ( Uj_jj/I^ ) sgn ( r ) , i— 1,2,3 (12) 
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If the state boundary conditions are such that a single 
principal axis (say, about axis 1 ) slewing is desired, then a 
solution satisfying the necessary conditions ( 10 ) and equations 
(5,7) for this problem can be obtained as 
qo(0)9^0, qiiOy^O, q2 ( 0 )=qa ( 0 )=0 , 
tJi(0)9^0, CJ2(0)= 6J3(0)=0; 

qi?^0, q2=q3“0^ 

Pl?^ 0 / P2=P3=0» 

COlj^O, C^2“ ^3=0; ri=0, r2=r3=0, 
ui=-(uib/Iii)sgn(ri) , U2=U3=0. 

The solution r 2 *r 3=0 implies a singular control problem because 
U 2 and U 3 can not be determined by Eq. (12). The possible 
existence of a singular solution in the general minimum time 
problem suggests that a unified method be needed to handle both 
singular and nonsingular cases. 

In some papers^'®, the integral of the sum of the squares 
of the torque components has been successfully used as a cost 
function, where there is no constraint on the control and the 
minimum time is not required. However, if some constraints on 
the control are added to this problem and the total slewing time 
is shortened sequentially, this problem may approach the minimum 
time control problem. These considerations motivate a successive 
approximation approach to solve the minimum time control problem. 
In this approach, an integral of a quadratic function of the 
controls is formally used as a cost function, i.e.. 


J=(l/2) 


tf 

/ u^Rudt 
0 


(13) 
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where R is a proper weighting matrix. It has been shown that^'®, 
for the case of rest-to-rest slewing with only 3 control inputs 
involved and with a longer slewing time used, the controls are 
approximately linear functions of time and do not reach their 
saturation levels. Therefore, when tf is shortened, some of the 
controls can be expected to reach their bounds and contribute 
more effort to the slewing. By the successive shortening of tf, 
a particular value, tf*, called the minimum time, can be 
obtained, during which either some (singular case) or all 
(nonsingular case) of the controls are of the bang-bang type. 

Apparently, in the approach described above, it is 
unnecessary to determine in advance whether the problem is 
singular or not and there is no need to determine the switching 
points as required in some other methods^®. 

Necessary Conditions 

The Hamiltonion for the system (5,7,13) is then 

H=(l/2)uTRu+pT5q+rT(i“l65lO +I"1 bu) ‘(14) 

where the costates, as before, satisfy the following necessary 
conditions for minimizing J, 

p=-(3HAq), or p=(l/2)wp (15) 

r=-(&H/3w), or r=g(w ,r )+(l/2 ) [qlp (16) 

where g( w , r) is a 3x1 vector function of co and r, and its 
detailed form can be found in Appendix A; [q] is a 3x4 matrix 
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-90 

-qo 

92 ■ 


[q] = 

q2 

93 

-qo 

“91 



.93 

-q2 

91 

-90 . 


The initial values of 

p and r 

are unknowns. 


The weighting matrix R in Eq. (13) is chosen as 

R»bTb (17) 

which is generally an nxn semi-positive-definite matrix, because 
the rank of the 3xn matrix, B, is assumed to be 3. Weighting 
matrices other than that given by Eq. (17) may also be possible 
candidates . 

Prom the necessary conditions (SH/du)=0, we have 
Ru+BTi“lr=0 
or 

B^BusB-B^I'lr (18) 

Premultiplying both sides of Eq. (18) by (BbT)“1b, one obtains, 

Bu*-I"lr 

By using the pseudo-inverse of the matrix B, B+, one can get u, 

u— B+I-lr»-BT(BBT)-l(i-lr) (19) 

The control laws are thenH 

'^i=“Uibsgn(B+I-lr)i, if |(B+I-lr )i| >uib; (20a) 

or 

Ui=-(B+I-lr)i, if |(B+I-lr)i| <uib (20b) 

1—1 ,2, ..., n. 

note that when B is a 3x3 nonsingular matrix, B"‘'=B“^. 
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A Linear Relation between q and p 

Before starting to solve the two-point boundary-value 
problem, it is useful to consider a relationship between q and 
p. It is already pointed out in Ref, 9 that 

p= p q 

where p is an arbitrary constant. However, one can find out that 
there does exist a linear relation between p and q in this 
problem, 

p(t)=Dq(t) (21) 

where D is a 4x4 constant matrix. To determine the constant 
elements of D, Eq, (21) is substituted into the differential 
equation for p in Eq.. (15), with the result, 

Dq=(l/2)CDq (22) 

The q in Eq. (22) is then replaced by Eq. (5), 

D(l/2)aq=(l/2)CDq (23) 

Since Eq. (23) is valid for arbitrary q, one has 

This relation is true for arbitrary values of <o only when D has 
the following form 



■do 

-di 

-d2 

-d3 

D* 

dl 

do 

-d3 

d2 


d2 

d3 

do 

-dl 


,d3 

-d2 

dl 

do 


and these constants, dg, di, d 2 / d 3 , can be determined by 
setting t=0 in Eq. (21), With the use of Eq. (24), the relation 
(21) can be rewritten as 
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Pol f<lO “<31 "^2 "< 33 ] [<^0 

Pi = qi <30 <33 "<32 <^1 

P2 <32 "<33 <30 <31 <^2 

P3j L<33 <32 "<31 qoj L<^3 

After substituting Eq. (25) into Eq. (16)/ there results 

r=g(cJ/ r)-(l/2)Cd (26) 

where d=[dx d2 dxl^z C is just the attitude matrix/ Eq. (1). It 
can be seen from Eq. (26) that r is independent of dQ. It is 
also true that u is independent of do because u depends only on 
r in Eqs. (20). This means that the arbitrary selection of do 
yields the same extremum control/ u. It is noted that a special 
choice of do can lead to the equivalent conditions considered in 
Refs. 3 and 9. With the use of Eq. (25)/ one can get 

do+<^l+<^2'^*^3“P^^® 1 (27) 

Since the choice of do is independent of the choice of d^/ 62 / 
and d 3 / a minimum value of the left side (hence the right side) 
of Eq. (27) is reached when do*0. This is the solution 
considered in Ref. 3. Also from Eq. (25) one can get 

do=pT(0)q(0) 

When do=0/ this equation gives a constraint on p(0). It is 
suggested in Ref. 9 that this constraint be used in the 
numerical iterations. But the choice of do other than zero is a 
more general result for this problem. It is not necessary to 
)ceep do*0 in each step of the computation. It is enough to keep 
one element of p(0) unchanged which is easier to use than the 
approach suggested in Ref. 9 , especially when q(0)j^[l 0 0 0]^. 
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Initial Values of Costates and the Slewing Time 

I 

1 Since the Euler rotation brings the attitude of the 

spacecraft from an initial quaternion to a final required 
quaternion through a simple rotation, it may take less time and 
consume less energy; it is reasonable to choose this rotation as 
a candidate for the starting solution of the iteration and hope 
that the optimal slewing is near the Euler rotation. This 
rotation will be called the "expected rotation", which is 

’ determined only by q( 0 ) and q(tf). 

The angular velocity and its derivatives for the Euler 
rotation can be expressed as 

co= 0 e, to » 0 e, * 0 e (28) 

where 0 (t) is the rotation angle and e=[ei e 2 63 is a unit 
vector along the the rotation axis (eigenaxis) which is 
determined by Eq. (4). Considering the analytical solution about 
a single principal axis maneuver in Ref. 3, 0 can be defined the 
same way about e, 

f 

} 6(t)*0(O)+0(O)t+(l/2)0(O)t2+(l/6)*0(O)t3 (29) 

► 

where 0 ( 0 ), 0 ( 0 ), 0 ( 0 ), and 0 ( 0 ) can be determined from the 

I 

boundary conditions of 0 (t) and 0 (t) at t= 0 , and t=tf. 

Without loss of generality, one can choose 

0(O)=O, 0(tf)=0* " (30a) 

where 0* is given in Eq. (4), and 

0(O)=0o, 0 (tf)= 0 f (30b) 
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The value of 0Q in Eq. (30b) needs to be determined from a given 
initial angular velocity vector, <o(0). Generally, this vector 
will not coincide with t3(0), the angular velocity of the Euler 
rotation at t=0, defined in Eq. (28), therefore, a difference 
vector, C, between them exists, 

£=0Qe-6)(O) 

Since only an approximate starting solution of the 

quasilinearization method is needed, it is enough to choose a Oq 

which minimizes By differentiating with respect to Gq 

and noting that e^e=l, one can get 

0O=eTco(O) (31) 

A similar derivation for 0f can be obtained. 

For the special case of 0f=O, the substitution of Eqs. (30) 
into Eq. (29) will result in 

0(O)=(60*/tf2)-(40o/tf) (32a) 

0*(O)=-(120*/tf3) + (60o/tf2) (32b) 

To approximately determine the initial values of p and r, 
Eqs. (7) and (26) are needed. By substituting u in Eq. (19) into 
Eq. (7) and solving for r, one can get 

r=IW icj -l2 <0 ( 33) 

r=(d/dt) (I OlfO )-i26J (34) 

At the same time, Eq. (26) can be rewritten as 

d=2cT[g( CO ,r)-r] (35) 
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Replacing to in Eqs. (33-35) by the relations (28-32) at t=0, 
one can get the approximate values of r(0) and d. p(0) can be 
determined by letting one of its elements equal a constant (say 
Po(0)=constant, if q(0)=[l 0 0 0]^) and by using Eq. (25) to 
solve for dQ and other elements of p(0). 

The starting solution needed in the quasilinearization 
algorithm may be obtained by integrating the differential 
equations (5,7,15-16,20) using the initial conditions p(0) and 
r(0) obtained above, as well as q(0) and (^(0). 

Initial Value of t<= 

Generally, to obtain the minimum time, one can always 
choose a longer slewing time, tf, at the beginning of the 
algorithm, and shorten it sequentially thereafter. But this may 
take more time, especially when how far the initial choice is 
from the real minimum time is not known. Therefore, a good 
initial value of tf being close to and larger than the minimum 
value is desired. For simplicity, only an estimation procedure 
of tf for the case in which B is a 3x3 unit matrix is discussed 
here. Suppose the slewing motion is an Euler rotation about a 
vector, e, through an angle, 0(t). By using the relations for 
in Eq. (28) into Eq. (6), one can get 

Ie0=02ele+u 

which can be expressed as the following 3 similar equations 

ai0=bi02+ciTi, i=l,2,3 (36) 

where af and bf are constants, Cf=ufi 3 , and -Cf is the normalized 
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control about the ith body axis and 

l*Ci|<l, i=l,2,3 

The 3 equations of Eq.(36) must be simultaneously valid for 
the same 6(t). Each of them with the boundary conditions, (30), 
can be considered as a minimum time control problem and solved 
analytically to obtained a minimum time (Appendix B). Since each 
of the minimum times for the associated equation means a lower 
bound of time during which the equation is solvable (no matter 
whether this equation is treated as a minimum time control 
problem or not ) , the largest one of these minimum times should 
be chosen as the initial value of tf used in the computation. 

A computation procedure has been developed which contains a 
series of cycles. The slewing time is chosen at the beginning of 
each cycle and fixed throughout the cycle. During each cycle, a 
quasilinearization algorithm called the method of particular 
solutions^2 ig used to solve the linearized state and costate 
equations. If this algorithm converges, a check is then made as 
to whether some (singular case) or all (nonsingular case) of the 
controls are of the bang-bang type. If yes, this slewing* time is 
designated the minimum time. If not, the assumed tf should be 
shortened and the next cycle begins. 

The numerical experience of using this procedure tells us 
that, for each cycle, the slewing time can not be made less than 
a certain value; in particular, it can not be made less than the 
real minimum time. Otherwise, the algorithm in each cycle will 
not converge. The closer the tf is to the real minimum time, the 
less shortening is required for the tf assumed in the previous 
cycle. 
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Numerical Results 


The methods described in the previous sections are applied 
to the SCOLE slewing motion. Pig. 1 shows the SCOLE 
configuration. It is composed of a Space Shuttle and a large 
reflecting antenna. The antenna is attached to the Shuttle by a 
flexible beam. Since only the motion of the rigid SCOLE is 
considered in this paper, the flexibility of the beam is 
ignored. The X, Y, and Z axes are the Shuttle axes corresponding 
to the roll, pitch, and yaw axes, respectively. The controls 
considered in this paper include three moments (ux“Ui,Uy»U 2 f 
U 2 =U 3 ) about the X, Y, and Z axes and two forces fy=U 5 ) 

applied at the center of the reflector in the X and Y 
directions. The inertia parameters of the SCOLE and the 
saturation levels of the controls are^: 

111=1132508, 122=7007447, 133=7113962, 

I12*“7555, 123=115202, Ii3=52293 slug-ft2; 

Uib=10000 ft-lb, i=l,2,3; 

Uib=800 lb., i=4,5 

The associated control influence matrix, B, in Eg. (6) is 
■ 1 0 0 0 130 ■ 

B= 0 1 0 -130 0 

.0 0 1 32.5 18.75. 

Some numerical results are presented in the following. 

(a) The singular case discussed in the previous sections 
for the SCOLE configuration without offset and for a symmetrical 
Shuttle Ii2=l23=0» For this case only the three Shuttle control 
torguers are used. The boundary conditions of the attitude are 
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such that the "expected rotation" is a single principal axis 
rotation. The computations show that when the slewing time 
approaches the minimum time, the control about the the slewing 
axis approaches the bang-bang type. The other two controls 
remain zero and there is no indication that they are going to 
make contributions to speed up the slewing. This result may 
imply that the singular solution is the time optimal solution 
for this case; otherwise additional control effort should 
participate in the slewing and a smaller slewing time should be 
obtained by using this algorithm. 

(b) The example of Ref. 3 is computed to test the method of 
determination of the initial costates. The results show that the 
guessed initial costates are very close to their converged 
values. To obtain the converged values (to seven digits), only 5 
iterations are needed in this computation. 

(c) The non-diagonal inertia matrix of the SCOLE and only 3 
controls (u^, U 2 r U 3 ) are used in this case. The expected 
rotation is a 20 degree, rest-to-rest rotation about one of the 
three spacecraft axes. Only the results for the "X-axis slewing" 
are given here because the results for the "Y-axis and Z-axis 
slewings" are similar. Pigs. 2 show the time histories of the 
controls and attitude angles (1-2-3 Euler angles) for this 
maneuver. Due to the non-zero offset of the inertia distribution 
of the SCOLE configuration, the controls, Uy and Ug, are no 
longer zero as they were in the singular case (a), but now are 
of the bang-bang type. The initial estimation of the minimum 
slewing time, tf, obtained using the method discussed in the 
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previous section , is tf ^ 0 )=12.5749 sec, which is very close to 
the minimum time, tf *=12.563034 sec, obtained in our computation. 

Pig. 3 shows the control torques for the "X-axis slewing" 
described in case (a) but for a slewing time, tf=15.37 sec, 
which is 2.8 seconds more than the minimum time, tf*. The 
controls are almost linear functions of time. ui is less than 
the saturation level; U 2 and U 3 are near zero. By comparing 
Pig. 3 with Pig. 2a, one can see that much more control effort 
(approximate 50%) is saved by using a little longer slewing 
time. Another feature of using a longer slewing time in the 
computation is that it needs fewer iterations for convergence 
than by using a shorter slewing time. These properties imply 
that, in the practical application of this problem, it is not 
necessary to seek exactly the minimum time, tf*, and the 
associated extremum controls. It may be enough to know the 
approximate values of tf* and the controls. 

(d) Pollowing case (c), two additional controls, U4 and U5, 
corresponding to the thrusters on the reflector are used. 

Pigs. 4 show the controls and attitude angles for the "X-axis 
slewing". The minimum time, tf *=3.9805382 sec, is greatly 
shortened as compared with case (c). Pigs. 5 show the controls 
and attitude angles for the "Z-axis slewing". The minimum time 
is tf *=15. 1441 sec. Unlike the case for the X-axis slewing, the 
attitude angle Gjj experiences a larger amplitude, though the 
expected rotation is about the Z-axis. This phenomenon is due to 
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the unsymmetric distributions of inertia about the X and Y axes. 
The closer the slewing is to the minimum time, the larger the 
amplitude of the 9x> 

(e) A general case is considered. Suppose the SCOLE is in 
an Earth orbit and the line of sight is to be directed toward 
the center of the Earth. The orbit coordinate system (x,y,z) is 
shown in Fig. 6. The initial attitude of the spacecraft is 
assumed as follows: the Y axis coincides with the orbital y 
axis, and the angular difference between the X and x (or Z and 
z) axes is «7 . 897224212 deg . The initial quaternion is, 
then, q( 0 ) = [cos ( «^/2 ) 0 sin(^/2) 0]^. According to Ref. 1, the 

unit vector along the line of sight in the rigid SCOLE 
coordinate system is 

RlOS*[ - 1112447 -.2410302 .9641209]T 
The direction cosines of the orbital z axis in the SCOLE system 
at the initial time are z/g=[sino^ 0 coso^jT. The angle between 
^LOS z/B at the initial time is GloS^O)* RloS*z/B= 20 deg . 
The eigenaxis of the expected rotation in the SCOLE system is 
determined by 

6“ ( Rlosxz/b ) / I Rlosxz/b 1 

The final required attitude quaternion can be obtained by using 
Eqs. (3-4). 

The guessed minimiam time for the case where only 3 controls 
are involved in this maneuver is tf =26. 3487 sec . This value is 
very close to the converged value, tf *=25. 003175 sec . It would 
be interesting if this estimation is compared with the result by 
solving the following classical minimum time control problem. 
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t 

I 


I 


I0-U, I'll^Umax 

where I is the moment of inertia about the principal line and u 
is the torque about that line, 0 is the angle of rotation. For 
the present case, the result obtained from the second method is 
tf=19.58 /sec . A possible explanation for the large difference 
here is that the classical problem greatly simplifies the 
inherent three dimensional nonlinear dynamics associated with 
the general SCOLE configuration. 

Figs. 7 show the controls and attitude angles for case (e) 
where U 4 and U 5 are also used. The minimum slewing time is 
obtained as tf *=8.691397 sec . The 0i,os Fig. 7b is the angle 
between the line of sight and the line of the target direction 
(from the spacecraft to the center of the Earth). 

Conclusions 

A useful numerical solution procedure for the minimum time 
attitude maneuver control problem of a rigid spacecraft has been 
developed and successfully applied to some practical examples. 

It can handle both the singular and nonsingular cases. It is 
shown through examples that the estimation methods used here for 
the initial costates and the minimtim slewing time are quite 
useful. The control profiles obtained in this paper may be 
useful for further research. 
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Appendix A The Term g( to, r) in Eq« (16) 

The term in Eq. (7) can be expressed as 

I-l col to=[PiG] w 

and 

= 2 2 2 

CO s[ (J-^ 0^2 ^3 ^2^3 ^3^ COj^6i)2]T 

where F and 6 are 3x3 matrices whose elements are constants 
associated with the inertia parameters of the spacecraft. Then 
the term rTl"^f3l«*) of the Hamiltonion, H, in Eq. (14) has the 
form 

h=rTi-l£5icj=[fi f2 f3 f4 fs feJT ^ 
where fi are 

[fl f2 f3l'^=FTr, [f4 fg felT^G^r 
The term g{ oi, r) in Eq. (16) can be obtained by 



2fi fg fs 


COl 

g( < 0 , r ) =“ Oh/Bto) =- 

f6 2f2 f4 


^2 

j 

fs f4 2f3 


0)3 


Appendix B Solution of Eg. (36) 

Eq. (36) can be rewritten as 

ai0=bi02+ciri (36) 

For simplicity, only the solutions for the following boundary 
conditions are considered here 

0(O)=O, 0(O)=O; 0(tf)=0*, 0(tf)*O (B-1) 

Suppose ai?^0, b^^^O and let b=bi/ai, c=ci/aj[ (suppose OO), one 
can rewrite Eq. (36) as 
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0=b02+cr 


(B-2) 


Since the control for this problem is of a bang-bang type with 
only one switching point f by integrating Eq. (B-2) and using Eq. 

( B-1 ) , one can get 

0=(c(e2be_i)/b]l/2, for r=l; (B-3) 

0-[c(l-e2b(0-0*) )/b]l/2, for.r=-l (B-4) 

By equating Eqs. (B-3) and (B-4), one can get 0=0s and 0=0s at 
the switching point , t=tg » 

0 s=(l/ 2 b)ln[ 2 /(l+e" 2 b 0 *) ] 

0s=tc(e2b0s-l)/bll/2 

Finally, by integrating Eqs. (B-3, 4) and using Eq. (B-1), one 
can get 

ts=cosh”^ ( e“^® s ) //-be , b<0; 
or 

tg=[ r^/2)-sin"^(e“'*^®s) l//bc, b>0 

and 

tf*tg+[ (■n/2)-sin“l(el5(0*-0s) )]//^, b<0; 
or 

tf*ts+cosh“l(e^^®*”®sM//Bc, b>0 
For the case 0(0 )^0, similar solutions can be obtained. 
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Fig. 1 Spacecraft Control Laboratory Experiment configuration 
( SCOLE ) 
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Fig. 2 X-axis slewing, tf =12. 563034 (s) 
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Pig. 3 X-axis slewing, tf=15.37 (s) 
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Pig. 4 X-axis slewing, tf=3. 9805382 (s) 


U 97 






(b) TIME (SEC) 

Pig. 5 Z-axis slewing, tf =15. 1441 (s) 
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SCOLE-example LOS slewing, tf=8. 691397 
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V. 


ASPECTS OF COMPUTATIONAL REQUIREMENTS FOR CONTROLLER IMPLEMENTATION 


V.I Discretization of Continuous Controllers: 


The full state variable feedback of the form 

U = KX (1) 

is usually proposed to control a linear plant given by 

X = AX + BU (2) 

using techniques such as pole placement^. Linear Quadratic Gaussian/ 

2 3 

Loop transfer recovery (LQG/LTR) , H techniques , etc. The imple- 

00 

mentation of the controller given in (1) and the generally needed 
state estimator are implemented using a microprocessor and thus 
the stability of the discretized control law as a function of the 
discretization step size is to be evaluated. 

To this end, the system djmamics and the control law as given 
by (2) and (1) are discretized using Euler's scheme as: 

X(i+1) = (I + A(A+BK)) X(i) (3) 

with 


U (i) = KX(i) 


(4) 


Balas^ has shown, using Lyapunov's stability criteria, that the stable 
continuous control law given by (2) is stable in the discrete time 
domain provided 


A < 


^min 

^max 



PA ) 
c 


(5) 


where 

Q = any positive definite matrix 
P is the solution of the Lyapunov's equation 
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PA + A + Q = 0 
c c 

X . ( ) , X ( ) are the minimxim and maximum eigen- 

min • max 

values of the matrices in the parenthesis 

A = A + BK 
c 

The disadvantage of the criteria (5) is its dependence on an 
arbitrary Q matrix and thus the absolute maximum value of A cannot 
be obtained. 


An improved stability criteria can be developed using the 
modal transformation on the closed loop system 



X = (A + BK)X 

(6) 

or 

• 

X = A X 
c 

(7) 

as 

X = T q 

(8) 


q = A q 

(9) 

where A = diag (X- , X«, ..., \) 
1 / n 



= ith eigenvalue of the matrix A^ and 
are assumed to be distinct. 

the eigenvalues 

The 

same transformation is applied to equation 

(3) as^ 


q(i+l) = t"^(I+ AA^)T q (i) 

(10) 

or 

* (1 + AXj^) 

(11) 


^ 1^2^...^ n. 

The £th difference equation is stable in the discrete time 
domain provided 



|1 + AX^l < 1 

(12) 

or 

(1 + AOj^)^ + (Aoij^)^ < 1 

(13) 
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where 


= o+ 2 tJi„ 


(14) 


i.e. 

A(20j^ + + oi^))< 0 


(15) 


thus, for stability, from the mathematical viewpoint, A either must 

2 2 

be negative and (2a + A(o + w )) must be positive which is 

X/ jL A/ 

practically unrealistic or, practically^ for positive A, 


A < 


-2a, 


2 2 


(16) 


It can be observed that A will be positive as a. is negative for a 
stable continuous control system. Thus the absolute maximum tolerable 
value of A is given by 


A 


= min ( 



) 


(17) 


on all £ where a^^= and o»^ is the undamped 

natural frequency of vibration. 

From the relation (17), it can be observed that the tolerable 
discretization step size. A, decreases with decreased damping ratio 
and increased maximum natural frequency of vibration as A can be 
approximated as (2^^/a)p . The stability of the continuous time con- 
trol system with repeated closed-loop eigenvalues is not considered 
here and will be attempted in the current grant period. 

Example: A typical large space structure^ with a maximum natural 

frequency of vibration, of IHZ (2Trrad/sec) is considered. The 
continuous controller is designed to provide 10% damping in all the 
modes and thus 
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( 18 ) 


Thus 


c. - 0.1 


\ ■ -'t 


W. - 


_ f 




A < 32 milli-seconds 
The closed* loop control system given by 
X = (A + BK)X 

can also be discretized as: 


(19) 


( 20 ) 


X(i+1) = X(i) (21) 

The eigenvalues of the matrix are the same as the 

eigenvalues of the similarity matrix T ^ e^^T where T and A are 
matrices given by the relations (8) and (9). As a result, the 
stability of the system (21) is decided by the eigenvalues of the 
matrix e^^ . Therefore, the magnitude of the eigenvalues of the 
matrix e^^ must be less than 1. 


i.e. 

1 

1 e ^ j < 1 

for all i 

or 

le'V 

< 1 for all 1 


And thus the following relation has to be valid for all 


( 22 ) 

(23) 


(24) 


The maximum A can be evaluated as the minimum of the A *s given in 

Xj 

relation (24) . The maximum A given by relation (17) is an approxi- 
mation of the value given by relation (23) . 
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V.II. Computational Requirements for Estimator and Controller; 

The estimator and the controller that are generally implemented using 
a microprocessor are described in the discrete time domain as^: 


X(i+1) 

= Ag X(i) + Y(i) 

(25) 

V(i) 

= K X (i) 

(26) 

Y(i) 

= CX(i) 

(27) 


where 

A 

X = n X 1 estimated state vector or a transformed 
state vector 

V = m X 1 input vector 

Y = £xl measurement vector 

A^,Be»K,C are appropriately dimensioned matrices. 

The minimum computational requirements are obtained through a modal trans- 
formation of the state vector, X, resulting in a diagonal system matrix A^. 

Thus the arithematic and analog to digital (A/D) conversion for sensors 
and digital to analog (D/A) conversion operations required to implement 
equations (17), (18) and (19) are as follows: 

No. of multiplications = n + nS. + nm 

No. of additions = n)t + nm (28) 

No. of A/D conversions = S. 

No. of d/a conversions = m 

Thus for a state of the art math-coprocessor such as Intel 80387 with 
25 MHZ clock and 10 KHZ sampling rate, the following are the typical times 
required for various computations. (Floating point arithematic is asstaned for 

g 

fastness and the range of values representable.) 
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Clock cycles 
(range) 


Average 

cycles 


Instruction 


Multiplication (C^^) 

32-57 

45 

FMUL 

Addition (C ) 
a 

29-37 

33 

FADD 

A/D conversion (A.) 

d 

0.1 msec 

- 


d/a conversion 

=0 

_ 

__ 


(negligible) 


Thus assuming that the number of actuators (m) and the number of 
sensors (£) are a fraction of the nximber of states (say £ = m = en) the 
computational time (C^) required to generate control signals from the 


continuous 

measurements is given by: 


thus 

C = n(l+2en) C + n(2en) C + en A^ 
t m ad 

(29) 


= n(l+2en) 1.8 + 2en^(1.32) + en(lOO) in y seconds 

(30) 

or 

= 6.24cn + n(1.8 + IOOg) 

(31) 


Thus the number of modes that can be handled by a microprocessor of the 
80387 t3^e to ensure stability of a large space structure with e = 0.25 
and A = 30 msec can be evaluated from the equation 

6.24 en^ + n(1.8 + IOOe) < 30x10^ ' (32) 


1.56a^ + 26. 8n - 30x10^ = 0 


-26.8 ±/( 26. 8) ^+4x1. 56x30x10^ 


2 X 1.56 


= 147.53 

The second root of the equation (33) being negative can be neglected. 
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The number of modes, being equal to half the number of states, that 
can be handled by a microprocessor of. the type 80387 is approximately equal 
to 74. Thus the exact discretization step size and the number of modes 
that can be handled by any given microprocessor can be arrived at by 
following the procedure given in this section and eliminating the assumptions made 
such as the actual times for A/D and D/A conversions, design safety factor 
for maximum et cetera. For a more stringent A. value, the relation 
given in (24) can be used. 
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VI. 


CONCLUSIONS AND RECOMMENDATIONS 


Mathematical models have been developed to predict the open and 
closed-loop d 3 mamics of the orbiting Spacecraft Control Laboratory 
Experiment (SCOLE), both during station keeping and, also during 
rapid slewing maneuvers of up to 20 deg. amplitude. It is seen that, 
in the presence of gravity-gradient torques, the system asstimes a 
new equilibrium position primarily due to the offset in the (flexible) 
mast attachment point to the reflector’s mass center. A robust 
control law based on an application of the linear regulator theory can 
be implemented for station keeping with maximum control efforts be- 
low saturation levels. 

For the slewing of a completely rigidized model of the SCOLE the 
LQR techniques can be extended to provide for relatively rapid slewing 
about each of the Shuttle's geometrical axes through amplitudes of up 
to 20 deg. These results can be compared with those provided by the 
numerical solution of the two point boundary value problem (TPBVP) 
associated with Pontryagin’s Maximum Principle, where the slightly 
faster slewing times are coupled with an increase in the over-all 
control effort. In connection with the latter approach, the minimum 
time attitude slewing of a rigid spacecraft has been examined with an 
integral of a quadratic function of the controls used as the cost 
function. Both singular and nonsingular problems can be treated in a 
unified manner. The resulting numerical solution to the TPBVP is 
based on a quasilinearization algorithm. General three dimensional 
slewing maneuvers (e.g. SCOLE antenna line of sight slewing) can be 
handled. 
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Extensions to the important problem of minimum time and near 
minimum time slewing will include the effects of flexibility in the 
subsequent grant year. A trade-off will need to be established between 
the rapidity of the slewing maneuver, control effort required by the 
variously placed actuators, and the ability to suppress the vibrational 
(flexible) amplitudes during and immediately following the slew. 

Finally, attention has also been focused on certain aspects of 
computational requirements for large space structural controller imple- 
mentation. An improved stability criterion is developed for a con- 
troller, designed to function in the continuous time domain, but which 
receives discretized observational inputs. An expression relating the 
maximum tolerable discretization step size to the undamped frequency 
and damping ratio of any mode in the continuous time system model has 
been developed. Future (on board) computational requirements are 
evaluated based on a current state-of-the-art microprocessor , assuming 
that the number of actuators and sensors are a selected fraction of 
the number of state components. Additional future attention to com- 
putational requirements should also be based on the mathematical models 
of the SCOLE system whose development has been completed during this 
grant period. 
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